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Abstract. We apply the method of comparison equations to study cosmological 
perturbations during inflation, obtaining the full power spectra of scalar and 
tensor perturbations to first and to second order in the slow-roll parameters. We 
compare our results with those derived by means of other methods, in particular 
the Green's function method and the improved WKB approximation, and find 
agreement for the slow-roll structure. The method of comparison equations, just 
as the improved WKB approximation, can however be applied to more general 
situations where the slow-roll approximation fails. 
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1. Introduction 

Anisotropies in the cosmic microwave background (CMB) radiation and inhomo- 
geneities in the large scale structures of the Universe have nowadays become a funda- 
mental tool to study the early universe [Q. Present and future data will allow us to 
discriminate among different inflationary models in the near future. For this reason, 
the comparison of observations with inflationary models requires theoretical advances 
in the predictions of the power spectrum of primordial perturbations beyond the low- 
est order in the slow-roll parameters e^'s, first obtained by Stewart and Lyth 2 (their 
definitions will be recalled in Section 

An analytic form for the full inflationary power spectra to second order in the 
slow-roll parameters was first obtained through the Green's function method (GFM 
henceforth) in Ref. Pj- In this way a characterization of the power spectrum to 
second order in the slow-roll parameters was given which was not just derived from the 
running of the spectral index (whose leading order is precisely 0{e^) j3j). An equivalent 
second order characterization has been recently obtained by means of the improved 
WKB approximation of Refs. [SJEl, which extended to second order in the slow-roll 
parameters previous results based on a more standard WKB approach |3|H]. The 
WKB approximation has confirmed the structure of the second order power spectra 
found with the GFM, within a numerical difference in the 0{e^) coefficients of 5% 
at most |5j. Compared with the GFM, the WKB approximation has the additional 
advantage that the slow-roll parameters do not have to be constant in time [HI and 
can therefore be applied to a wider class of inflationary models. 

The purpose of this paper is to illustrate the use of the method of comparison 
equations (MCE in brief) [HI EH ^3 to predict inflationary power spectra. We shall 
see that this method yields exact results for the case of constant slow-roll parameters 
(e.g. power-law inflation) and polynomial structures to second order in the e^'s in 
agreement with the GFM and WKB approximation. The MCE however has the 
advantage of being more accurate to lowest (leading) order, whereas other methods 
(namely, the WKB 0, improved WKB and GFM 0) reach a similar accuracy 
to next-to-leading order. We shall also discuss cases for which our present method 
appears more flexible than the slow-roll approximation. 

In the next Section we briefly review the general MCE and in Section|31the theory 
of cosmological perturbations. We then apply the MCE to cosmological perturbations 
in Section 0] where we also analyze the error around the "turning point" in detail. 
In Section El we analyze power-law inflation, chaotic inflation and the arctan model; 
we expand our general results to second slow-roll order and compare with analogous 
results obtained with the GFM. We finally comment on our results in Sectional Some 
more technical details are given in two Appendices. 

2. Method of comparison equations 

Let us briefly review the MCE (the name is due to Dingle ^|). It was independently 
introduced in Refs. [UJ llOj and applied to wave mechanics by Berry and Mount in 
Ref. The standard WKB approximation and its improvement by Langer 12 are 
just particular cases of this method and, recently, its connection with the Ermakov- 
Pinney equation was also studied 
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Let us consider the second-order differential equation 



(1) 



where is a (not necessarily positive) "potential" (or "frequency" ) , and suppose that 
we know an exact solution to a similar second-order differential equation, 



d(T2 



U{a) = 



(2) 



where Q is the "comparison function" . One can then represent an exact solution of 
Eq. in the form 

\ -1/2 

Uia) , (3) 



dx ) 



provided the variables x and a are related by 



tj2(a;): 



dcrV 
d^ / 



da\^/' d^ 



da; / Ax^ \Ax 



dcr 



-1/2 



(4) 



Eq. can be solved by using some iterative scheme, in general cases ^1 El or for 
specific problems jl5lll()| . If we choose the comparison function sufRciently similar to 
w, the second term in the right hand side (r.h.s.) of Eq. Q will be negligible with 
respect to the first one, so that 



(5) 



On selecting a pair of values x^ and cto such that ctq = the function <j{x) can 

be implicitly expressed as 



J Xn J (Tn 



(6) 



where the signs are chosen conveniently |. The result in Eq. (|2Jl leads to a uniform 
approximation for xi^): valid in the whole range of the variable x, including "turning 
points" §. The similarity between Q and lu is clearly very important in order to 
implement this method. 



3. Cosmological perturbations 



Let us begin by recalling that scalar (density) and tensor (gravitational wave) 
fluctuations on a flat Robertson- Walker background with scale factor a 

dr^+r^dn^) , 



ds^ = {-dr]^ 



(7) 

are given respectively hy fi — fig = aQ and fi — fiT = ah, where Q is the Mukhanov 
variable |17[ [T5| and h the amplitude of the two polarizations of gravitational 
waves ^1 120]. The functions /z must satisfy the one-dimensional Schrodinger-like 
equation 

j2 



d- 

drj-^ 



/I = 



(8) 



X We recall that ^{x) is the same quantity as used in Ref. ISIISIIHI . 

§ With this term, borrowed from point particle mechanics, one usually means a real zero of the 
"frequency" uj. 
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together with the initial condition (corresponding to a Bunch-Davies vacuum) 



-ikrj 



/2k 



(9) 



In the above rj is the conformal time (derivatives with respect to it will be denoted by 
primes), k the wave-number, H = o! jc? the Hubble parameter and 



VL\k,r]) EE k^ 



(10) 



where z = z-g, = a} 4>' / H for scalar and z = z^ = a for tensor perturbations (0 
is the homogenous inflaton). The dimensionless power spectra of scalar and tensor 
fluctuations are then given by 



fc3 

2 7r2 



and the spectral indices and runnings by 



_ dlnP^ 
^ ^ " dlnfc 



Tlx 



k=k. 



d^ \nV, 



as 



dlnfc 
d^ Inn 



(dlnfc)2 ^^^^ ' ""^ - (dlnA:)2 

where fc* is an arbitrary pivot scale. We also define the tensor-to-scalar ratio at fc = fc* 
as 



(11a) 



(lf6) 



(11c) 



k=k. 



(lie 



k=k. 



Finally, in the following wc shall often make use of the hierarchy of horizon flow 
functions (HFF in short, also referred to as slow-roll parameters) e^'s defined by 



ei 



H 



Cn+l = 



He,, 



n>l 



(12) 



where dots denote derivatives with respect to the cosmic time clt = a drj. 



4. MCE and cosmological perturbations 

In order to apply the MCE to cosmological perturbations, we shall start from the 
same equation as was used with the improved WKB method in Refs. |3IH1, to which 
we refer for more details. We recall here that the WKB approximation can be more 
effectively applied after the following redefinitions of the wave-function and variable, 

X = {l-,i)"'<:-'"li (13a) 

. = 1.(J-). ,13,, 

This yields an equation of the form ^ with the "frequency" n{k, rj) of Eq. ((H)) replaced 

by 

[1 - ei{x)\ 
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with given, respectively for scalar and tensor perturbations, by 



1 

4 



3-£i 
3-£i 



- (3-2£i)e2 (l-2£i)£2e3 (l-4£i)e^ 
2(l-ei)2 + 2(l-ei)3 + 4(l-ei)4 

2 _ _ _ _ _ /r> , _ N _ _2 



£1 £2 



£162 £3 (2 + ei)eie2 



(15a) 



(156) 



2(l-ei)2 2(l-ei)3 4(l-£i)4 ' 

where we omit the dependence on x in the for the sake of brevity. The point x = Xq 
where the frequency vanishes, oj{xo) = (i.e. the classical "turning point"), is given 
by the expression 

Xo=ln[P (l-£i)] , (16) 

where we have defined D = i'Ixq) and £i = £i(a:o)- We now choose the comparison 
function 

„2 a 



(1-f-l 



(17) 



and note that then ctq — xq. Solutions to Eq. can now be expressed, by means of 
Eqs. lO, © and as 



/e(a) 

uj{x) 



J±L> 



1-ei 



(18) 



where the J's are Bessel functions ^21 j and the initial condition ^ can be satisfied by 
taking a linear combination of them. However, in contrast with the WKB method [71|S] 
and as we pointed out in Section[21 MCE solutions need not be matched at the turning 
point, since the functions (|18|) are valid solutions for the whole range of the variable 
X. Eq. © at the end of inflation, x — Xi, becomes 



^Xi) ~ -e(af) 



In 



v - 9(crf) 



1 + ln 



1 -ei 



- In (2 v) 



(19) 



where the super-horizon limit Xf <C 2^0 (ff 
One then has 



2v 

— exp 



C(a^f) 



-oo) has been taken in the second line. 



(20) 



Finally, on using Eq. lf2n|) , we obtain the general expressions for the power spectra to 
leading MCE order, 

3 



TT £1 TOpj \aH J (1 — £1) 



.2 ?s 



16i?2 



'PI 



aH J (1 — £1) cjt 



(21a) 



(216) 



where mpi is the Planck mass and the quantities inside the square brackets are 
evaluated in the super- horizon limit, whereas the functions 



[1 - cos (2 TT I?)] [r(l 



(21c) 
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describe corrections that just depend on quantities evaluated at the turning point 
and represent the main result of the MCE applied to cosmological perturbations |j. 
The expression in Eq. (|21c|) is obtained by simply making use of the approximate 
solutions H18|) and their asymptotic expansion at a; — s- oo to impose the initial 
conditions In the WKB calculations one finds a similar factor but, in that 

case, using the Bessel functions of order 1/3 leads to a large error in the amplitudes. 
The MCE instead uses Bessel functions of order D, with P = 3/2 to leading order 
in the HFF (i.e. the right index for the de Sitter model), which yields a significantly 
better value for the amplitudes of inflationary spectra. 

The MCE allows one to compute approximate perturbation modes with errors in 
the asymptotic regions (i.e. in the sub- and super-horizon limits) which are comparable 
with those of the standard (or improved) WKB approximation |3 |H]. Since these 
methods usually give large errors at the turning point 8 (which produce equally 
large errors in the amplitude of the power spectra) it will suffice to estimate the error 
at the turning point in order to show that the MCE is indeed an improvement. To 
leading order (that is, on using the approximate solution (|18|l 'l. the MCE gives an 
error at the turning point of the second order in the HFF, which means that we have 
a small error in the amplitudes of the power spectra. Unfortunately, this error remains 
of second order in the HFF also for next-to-leading order in the MCE. We shall see 
this by applying Dingle's analysis ^U] for linear frequencies to our case (|14|) . We start 
by rewriting Eq. as 



(l-e-i)2 



3a| 



1 Ojj 

2 (71 







(22) 



where we dropped the x dependence in to and a and the order of the derivatives is 
given by their subscripts (cti = dcr/dx, uif = dcu'^/dx, etc.). Note that the term in 
square brackets is the error obtained on using the solutions H18() . We then evaluate 
Eq. and its subsequent derivatives at the turning point (i.e. a.t x — xq), 



{uj'-ale''{a)} + 



3al 
4 a? 



1 CT3 

2 (71 







(23a) 



2(72 <yi (a) 



2e 



2 a ^3 



(73 (72 3 crS 



1 (74 

2 a\ 



{23b) 



2 ((73(71 



2 af 



^1) 



9 crj 17 as CT2 



((7) - 

5 (74 (72 

" 2 



2(7f e 



2 _2 cr 



(l-6l)^ 



(2(7^+5(72) 



1 CTs 

2 ^ 



= 



(23c) 



2 ((74 (71 + 3(72(73) ((7) 



2(7 



87 (73 erf 

y 



2 (7i e' 
(l-^-i) 

2 , 



(4(7^ + 18(72 (7^+7(73 (71 + 12(7^ 



18 



4 _ 27 (74 (7| _ (7| (72 g5 ^2 



13 (73 (74 
Y 



1 CT6 

2 (71 



0(?3(f) 



'1 '^l ^ "1 ^1 "1 ^ "1 

II Inside the square brackets we recognize the general results given by the WKB approximation 171 151 . 
The "correction" g{xo) accounts for the fact that in Eq. 1171 is a better approximation than 
Langer's 151 IT^. 
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where 0^ (a) was defined in Eq. H17|) and we omit two equations for brevity. In order 
to evaluate the error at the turning point 



Atp — 



4 a? 



1 03 

2 ^ 



(24) 



we ignore the terms in square brackets and equate to zero the expressions in the curly 



brackets in Eqs 
cr = In [(1 — ei) P 

1/3 



Q3m and so on. This leads to 



0-2 = 



^1 
2P2 

1 



5 (2P2) 



1/3 



6 (21/3 



ox2 



2d 



2/3' 



<73 = 



3-21/3^^2 



16 ujj 



175 (p2/5 ^2)5/3 25 (i^4^2)l/3 175 p2 7(31/2^^2) 



2/3 



(25 a) 
(25&) 

(25c) 
(25 rf) 



and similar expressions for (74, (75 , and erg which we again omit for brevity. On inserting 
Eqs. H14() . H15a|l and ()155|l in the above expressions, we find the errors to leading MCE 
order 



32 

315 
32 

315 



£1 £2 



£1 £2 



22 
315 



£2 £3 



(26 a) 
(26 &) 



for scalar and tensor modes respectively. 

On iterating this procedure we can further obtain the errors for the next-to-leading 
MCE order Ai^^p. We first compute next-to-leading solutions to Eqs. H23af) - (l23dl) and 
so on by inserting the solutions found to leading order for tri , cr2 , . . . , ce into the 
corrections (i.e. the square brackets) and into all terms containing PO]- This 

leads to 



^(1) 

^(1) 
^TP,T 



31712 

331695 
31712 



21598 
331695 



£1 £2 - £2 £3 (27a) 

331695 ^^^^' 

which show that the next-to-leading MCE solutions lead to an error of second order in 
the HFF, too. We suspect that this remains true for higher MCE orders, since there 
is no a priori relation between the MCE and the slow-roll expansions. Let us however 
point out that the above expressions were obtained without performing a slow-roll 
expansion and therefore do not require that the be small. 



5. Applications 

In this section we apply the formalism developed in the previous section to some 
models of inflation. We shall expand our general expressions H21a|l - (|21cf) to second 
order in the HFF and compare them with other approximation methods used in the 
literature. 
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5.1. Power-law inflation 

In this model \2'6l |^ , the scale factor is given in conformal time by 

a{v)^io\v\'-"' , (28) 
where (3 < —2 and io = corresponds to the (constant) Hubble radius for de Sitter 
(/3 — —2). Since the HFF are constant, 

2 + /3 



71 > 1 



1 + ' 

the MCE yields the exact power spectra, spectral indices and runnings. 



Vh 



16 £■ 



PI 



./(/3)fc^^+ 



where £pi = m^l is the Planck length and 

1 



i.pi 

TT 



22/3+1 [i_cos(2 7r|/3- 



r2 {13 ■ 



1 

TT 



r(/3 + 



H 2 



2/3+1 



(29) 

(30 a) 
(306) 

(30c) 



with r the Gamma function. The spectral indices are rig — 1 = nx = 2/3 + 4 and their 
runnings ag = a-r = 0. Finally, the tensor-to-scalar ratio becomes 
2 + 13 



i? = 16 



1 + 13 ' 



(31) 



which is constant as well. 



5.2. Leading MCE and second slow-roll order 

We now consider the results H21all - (|21c|) given by the MCE to leading order (denoted 
by the subscript MCE) and evaluate them to second order in the HFF (labelled by 
the superscript (2)) for a general inflationary scale factor. A crucial point in our 
method is the computation of the function ^ defined in Eq. © which can be found in 
detail in Section HI of Ref . [HI . For the sake of brevity, we shall not reproduce that 
analysis here but it is important to stress that, in contrast with the GFM and other 
slow-roll approximations, it does not require a priori any expansion in the HFF since 
Eq. (34) of Ref. [H] is exact and higher order terms are discarded a fortiori. From that 
expression, upon neglecting terms of order higher than two in the HFF, we obtain the 
power spectra 

'PfL.=-^U-nC+l)e.~Ce,+ (2C' + 2C+^^^^\ el+ (]-C^ + ^^ -1 



2C2-2CZ?MCB+-DMCB-C-2Cln(2) + 2i^MCB ln(2) + ^-y^ ei £2 
-C Z?MCB + \ Z^Lb - C ln(2) +Dmcb ln(2) + ^ - £2 £3 
^ [- 2 £1 - £2 + 2 (2 C + 1 ) £2 + (4 C - 2 i?MCB - 1 ) £1 £2 + C £2 - Z?MCB £2 £3] In ^ 
-i(4£2 + 2£i£2 + £^-£2£3)ln2('AU (32a) 
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-2Ci5MCE + i^MCE ~2C~2C ln(2) + 2i?MCB ln(2) + — - — ) £162 



12 9 



[-2ei + 2 (2C + 1) £2 _ 2 (Z^mce + 1) eie2] In 



k 



k 

].[Ael-2e,e,)\n'{^]] , (326) 



where I^mce = - ln3 « -0.7652 and C = In 2 + 7e - 2 « -0.7296, with 7e the 
Euler-Mascheroni constant. A clarification concerning the g{xoys is in order. Since 
the turning point does not coincide with the horizon crossing where the spectra are 
evaluated [S], we have used the relation 

€i{xo) ~ei-ei e^+i In , (33) 

in order to express the g(xo)'s as functions of the crossing time (HFF with no explicit 
argument are evaluated at the horizon crossing). The spectral indices Ijll6|l are then 
given by 

4'LcE - 1 = -2 ei - £2 - 2 - (2 Z^mce + 3) ei £2 - ^mce £2 £3 (34a) 
^T.lcE = -2 £1 - 2 e? - 2 (Dmce + 1) £1 £2 , (346) 



and their runnings Hllc|l by 

4^LcE = -2£i£2 - £2£3 (34c) 

«T,lcE = -2eie2 . (34rf) 
The tensor-to-scalar ratio (llldjl becomes 

i?^^^ / 7r2 \ /I 



= 1 + C£2+ — + 5 £i£2+ - 



16ei V 2 / " " V2 



^2 



C i^MCE - I Dl,^ + C ln(2) - i?MCE ln(2) - ^ + 1^ £3 . (35) 

The polynomial structure in the HFF of the results agrees with that given by the 
GFM [3 121] and the WKB approximation |S1 E] (the same polynomial structure 
is also found for the spectral indices by means of the uniform approximation j25j). 
Let us also note other aspects of our results: first of all, the factors 5(^0) modify 
the standard WKB leading order amplitudes [3 |H] so as to reproduce the standard 
first order slow-roll results; secondly, we found that C and I^mce are "mixed" in the 
numerical factors in front of second order terms (we recall that Dmce differs from C by 
about 5%); further, Dmce — D^jkb of Refs. |3IH1E1E1- The runnings as and ax are 
predicted to be ©(e^) Jl], and in agreement with those obtained by the GFM p^l26|. 



5.3. Second slow-roll order MCE and WKB versus GFM 

We shall now compare the slow-roll results obtained from the MCE and other methods 
of approximation previously employed H] with those obtained using the GFM in the 
slow-roll expansion. We use the GFM just as a reference for the purpose of comparing 
the other methods among each other and of illustrating deviations from pure slow-roll 
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Figure 1. Percentage differences I36i between scalar-to-tensor ratios given by 
the GFM and those obtained from the WKB, WKB* and MCE. 



results. For a given inflationary "observable" Y evaluated with the method X, 
denote the percentage difference with respect to its value given by the GFM as 



we 



100 



Y - Y 
Y 



% 



(36) 
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I^Smce - "ScfmI X 100 AP^MCE(fc*) 



Figure 2. Left box: absolute difference between scalar spectral indeces ng 
evaluated with the MCE and GFM (rescaled by a factor of 100 for convenience). 
Note that relevant differences of order 0.05 (shown as 5) occur at the boundaries 
of the intervals considered for the e^'s. Right box: percentage difference I36i for 
the amplitudes of scalar perturbations at the pivot scale k = k,t between the 
MCE and GFM. See Fig. 0for the meaning of dot size. 




Figure 3. Percentage difference I36i between the MCE and the GFM for 
P^(fc*) restricted on the hypersurface £2 = — 2ei, which corresponds to a scale- 
invariant spectrum to first order in the slow-roll expansion. The graph is given 
for < ei < 0.5 and -0.5 < £3 < 0.5. 

where X = WKB stands for the first WKB and second slow-roU orders X = WKB* 
for the second WKB and second slow-roll orders |S] and, of course X = MCE for the 
result obtained in this paper. Note that for the case X = WKB* we shall set the three 
undetermined parameters 5s = 6t = c^s = 2 in order to minimize the difference with 
respect to the results of the GFM (see [Appendix A| and Ref. (HI)- 

In Fig.^we show, with dots of variable size, the percentage differences H36|l for the 
scalar-to-tensor ratios R at the pivot scale fc = fc* for < ei < 0.5, |e2| and jeaj < 0.5. 
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From the plots it appears that the level of accuracy of the MCE is comparable to 
that of the WKB* and both are (almost everywhere) more accurate than the WKB. 
However, the MCE achieves such a precision at leading order and is thus significantly 
more effective than the WKB*. In Fig.[21we show the difference in ns and the relative 
difference in P(j(fc*) between the MCE and GEM. In Fig.|2|we finally plot the relative 
difference in Pi; (A:*) for £2 = — 2ei (the scale invariant case to first order in slow- roll 
parameters) . 



5.4- GFM and slow-roll approximation 

Let us further compare our results with those obtained by the GFM. As we stated 
before, the most general result of the present work is given by the expressions for the 
power spectra H21q|) and (|21&|l with the corrections (|21c|l . In fact, the spectral indices 
and their runnings are the same as found with the standard WKB leading order [3 EI ■ 
On the other hand, the main results of the GFM are the expressions for power spectra, 
spectral indices and runnings ". . . in the slow-roll expansion." that is, for small HEE 
(see, e.g., Eqs. (41) and (43) in Ref. 3 ). Here, and in some previous work 171, 15!, we 
instead obtain results which hold for a general inflationary context, independently of 
the "slow- roll conditions" ^ or "slow-roll approximation" +. The slow- roll case is just 
a possible application of our general expressions which can be evaluated for any model 
by simply specifying the scale factor. Our general, and non-local, expressions in fact 
take into account all the "history" of the HEE during inflation. For example, we note 
that the MCE at leading order reproduces the exact result for power-law inflation (see 
Sec. I5.1|l . whereas the GEM reproduces those to the next-to-leading order 3 . 

Chaotic inflation is another example for which a Taylor approximation of the 
HEE (such as the one required by the GFM 3 or our Green's function perturbative 
expansion in [Appendix B| ) may lead to inaccurate results. We consider a quadratic 
potential 

y(0) = im2 02^ (37) 

where m is the mass of the inflaton 0. In a spatially flat Robertson- Walker background, 
the potential energy dominates during the slow-rollover and the Friedman equation 
becomes 



3 m|,j 



2 i2 

m 



3 m|,j 



(38) 



The Hubble parameter evolves as 

H = ~^ll? . (39) 

On using the equation of motion for the scalar field 

+ 3i/0 + m20 = , (40) 

and neglecting the second derivative with respect to the cosmic time, we have 
3 7J ~ - 0. Eq. (jSHl then yields 

2 

H^-^^Ho, (41) 

1 Our approximation method does not require small HFF. 
+ We expand instead according to the scheme given in Ref. 



Method of comparison equations for cosmological perturbations 



13 




Figure 4. Numerical evolution of H (solid line) and its analytic approximation 
(dashed line). The initial condition corresponds to Hi ^ 8.2m (time in units of 
1/m) when the inflaton starts in the slow-roll regime with a value 4mpi. 



which leads to a hnearly decreasing Hubble parameter and, correspondingly, to an 
evolution for the scale factor which is not exponentially linear in time, i.e. 

H(t) :^ Ho + Hot , (42a) 

a{t) ~ a(0) exp (^Ho t + Ho*—^= a(0) exp [N{t)] . (42&) 

In Fig. ^ the analytic approximation (|42all is shown to be very good on comparing 
with the exact (numerical) evolution. We are interested in the HFF in the slow-roll 
regime. From the definitions (|12|l we then have 

^lit)^-- ^^^^ (43a) 

Hot 



{Ho 

(^Ho + Hot) 



^n{t) = -- . ^2ei(0 , n>2, (43&) 



which we plot in Fig. |S1 

Let us take the end of inflation at fond — 24/m (so that the analytic 
approximation (|42q|) . (|42&|) is very good till the end), corresponding to iVond = 125 e- 
folds. We then consider the modes that leave the horizon 60 e-folds before the end of 
inflation and take that time (i = ~ 8/m or N^, = 125 — 60 = 65) as the starting 
point for a Taylor expansion to approximate the HFF * , 

e, (AiV) ~e, + e, e,+i AiV + ^ e, {ej^^ + e,+i e.+s) AiV^ , (44) 

where the e^'s in the r.h.s. are all evaluated at the time and AA^ = N — N^,. The 
percentage error with respect to the analytic expressions H43a|l and (|436|) , 

. 100 % , (45) 

* One can of course conceive diverse expansions, e.g. using other variables instead of A'^, however the 
conclusion would not change as long as equivalent orders are compared. 
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Figure 5. HFF in the slow-roll regime: the dashed line is ei and the solid line is 
all the e„'s with n>2 (time in units of l/m). 




Figure 6. Percentage error for ei to first and second order in AN (dashed and 
solid line respectively; time in units of l/m). 



with AN — N{t) — N^, is then plotted in Fig. Elfor i = 1. It obviously becomes large 
rather quickly away from t — t^ and we can immediately conclude that, had we used 
a Taylor expansion to approximate the HFF over the whole range of chaotic inflation 
(—65 < A7V < 60), we would have obtained large errors from the regions both near 
the beginning and the end of inflation. 

In general, we expect that the Taylor expansion of the HFF will lead to a poor 
determination of the numerical coefficients (depending on C in Eqs. H32Qjl and (|32&(l 'l in 
the second slow-roll order terms for those models in which the HFF vary significantly 
in time. In fact, we can calculate the integral of Eq. (21) in Ref. |3] with ei instead of 
the complete g (In a;) and yo{u) instead of y{u). The percentage difference between this 
integral calculated using the Taylor expansion (|44|1 with respect to the same integral 
calculated with ei in Eqs. H43 a|l and l|43 6|l is 92 % and 88 % to first and second orders 
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Figure 7. Spectrum of the Mukhanov variable Q (Pq = fc^ |Qfe p/(2 tt^)) 
evaluated at the end of inflation for the chaotic model of Eq. I37i . The dots 
represent the numerical values and the solid line the analytic fit based on 
Eqs. 134 al and 134 cl obtained by the second slow-roll order MCE approximation. 
The first order and GFM analytic results are not shown since they are almost 
indistinguishable from the MCE result plotted, fc* crosses the Hubble radius at 
</>* ~ Srripi (the lines are normalized at Vfi-^^ fc*). 



in AiV, respectively. The MCE method does not require such an expansion and does 
therefore not suffer from this restriction jj. 

We have also compared our analytic MCE results in Eqs. (|34Qfl and (|34cl) with 
the numerical evaluation of the spectrum for scalar perturbations. In particular, we 
have evolved in cosmic time the modulus of the Mukhanov variable Q (recalling that 
Pq = fc"^ |Qfep/(2 TT^)), which satisfies the associated non-linear Pinney equation j28j . 
The initial conditions for the numerical evolution are fixed for wave-lengths well 
within the Hubble radius and correspond to the adiabatic vacuum, i.e. \Qk(ti)\ = 
l/{a{t,)Vk) and \Qk\iti) = -H{t,) \Qk{ti)\ The agreement of the analytic MCE 
approximation with the numerical results is good, as shown in Fig. [T] 



5.5. Beyond slow-roll 

The slow-roll approximation is quite accurate for a wide class of potentials. However, 
violations of the slow-roll approximation may occur during inflation, leading to 
interesting observational effects in the power spectra of cosmological perturbations. 
An archetypical model to study such violations is given by the potential 



2 

1 arctan N ■ 



(46) 



introduced with iV = 5 in Ref. 30 . In Fig. |51we show that the MCE result expanded 
to second slow-roll order provides a very good fit for the power spectrum of scalar 
perturbations even in situations where the slow-roll parameters are not very small 
(see Fig. ISJ . This example also shows that second slow-roll order results are much 
better than first order ones (analogous results were also obtained with the GFM in 
Ref. |2ni)- 

(t We again refer the reader to Section III of Ref. ^ for more details. 
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Figure 8. Evolution of with the value of the inflaton (in units of mpi) in the 
arctan model of Eq. I46i : ei (solid line), £2 (long-dashed line) and 62 /H = £2^3 
(short-dashed line). 
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Figure 9. Spectrum of the Mukhanov variable Q (Pq = [Qj. p/(2 7r^)) for the 
arctan model of Eq. 1461 . The dots represent the numerical values, the solid line 
the analytic results from Eqs. 134 ai and 134 ci obtained by the second slow-roll 
order MCE approximation and the dashed line those given by the first order slow- 
roll approximation. fc« crosses the Hubble radius at 0, ~ — O.Smpi (the lines are 
normalized at lO'^^^^k,) and the spectrum is evaluated at ~ 55 e-folds afterwards. 



6. Conclusions 

We have presented the apphcation of the method of comparison equations to 
cosmological perturbations during inflation. By construction (i.e. by choosing 
a suitable comparison function), this approach leads to the exact solutions for 
inflationary models with constant horizon flow functions e^'s (e.g. power-law inflation). 

The main result is that, on using this approach to leading order, we were able 
to obtain the general expressions H21all - (|21c|l for the inflationary power spectra which 
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are more accurate than those that any other method in the hterature can produce at 
the corresponding order. In fact, the MCE leads to the correct asymptotic behaviours 
(in contrast with the standard slow-roU approximation pi) and solves the difBculties 
in finding the amplitudes which were encountered with the WKB method [Tj. 

Starting from the general results H21q|) - (|21c() . we have also computed the full 
analytic expressions for the inflationary power spectra to second slow-roll order in 
Eqs. H32a|l and (|32bll and found that the dependence on the horizon flow functions 
Ci's is in agreement with that obtained by different schemes of approximation, such as 
the GFM 131 and the WKB approximation Moreover, the results obtained with 

the MCE do not contain undetermined coefficients, in contrast with second slow-roll 
order results obtained with the WKB* 0. 

Let us conclude by remarking that, just like the WKB approach, the MCE 
does not require any particular constraints on the functions e^'s and therefore has 
a wider range of applicability than any method which assumes them to be small. 
As an example, we have discussed in some detail the accuracy of the MCE for the 
massive chaotic and arctan inflationary models. We have shown that the MCE leads 
to accurate predictions even for a model which violates the slow-roll approximation 
during inflation. 
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Appendix A. Impact of undetermined parameters in the WKB* 
approximation 

To clarify the impact of the undetermined parameters bg, b^, and ds = 2 on the 
results to second order in the HFF, the percentage differences between the numerical 
coefficients of some second order terms in the power spectra are displayed in Table IXT1 
and of some of those in the scalar-to-tensor ratio in Table IA2I All the entries are for 
bs — b^ — dg ~ 2, as is used in the main text, since other choices would give larger 
differences, as we explicitly show in Fig. lAll for the numerical coefficient in front of 
ei 62 in Ph (since the behaviour is the same for all coefficients, we do not display other 
plots). 





e 'tPc Ph 


4-:P( 


eie2; Pcfe =2) 


ei £2; Ph{bT = 2) 


WKB 


2.3% 


3.0% 


353.3% 


37.2% 


WKB* 


0.03% 


0.06% 


29.1% 


2.0% 



Table Al. Percentage differences between the numerical coefficients of some 
second order terms in tlie power spectra given by WKB and WKB* with respect 
to those obtained by the GFM. 



Appendix B. Perturbative approximations 

In this Section we present a possible method of obtaining the complete result to second 
order in the HFF. We start by noting that any given function / = f{x) = /(e, x) can 
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£2 £3; {ds = 2) 


ei £2; {bs - 6t = 0) 


4 


WKB 


56.7% 


1.5% 


117.3% 


WKB* 


1.3% 


0.03% 


1.5% 



Table A2. Percentage differences between the numerical coefficients of some 
second order terms in the scalar-to-tensor ratio given by WKB and WKB* with 
respect to those obtained from GFM. Note that the numerical coefficient of €2 ea 
in R is the same as that in P^. Further, the coefficient of ei £2 does not depend 
on any undetermined parameters if 63 = 6t- 



be written as the power series 



n=0 



Tl = 



where = ei(x») are the HFF evaluated at the horizon crossing (i.e. at x — x^) and 
f{ei, x) — f{€i — ii,x). On using Eqs. HB.1|) to express both uj"^ and x in Eq. we 
obtain 



Xo{x) V —ru;l{x) + 
, n! 

n— 1 

where we have used 

2/ 



+ — r^mW 



m=0 



m! 



Xo(a;) = 



™7t 



(B.3) 



with (^q{x) — uj^iti — ii',x) and Xoi^) is a hnear combination (i.e. the MCE leading 
order solution) of Eq. p8|l for <d^{x) ~ (^o{x) and cr ^ a; in the argument of the Bessel 
functions. One may solve Eq. I)B.2|I for higher orders by introducing 

a;" 

— r Xnix) = g„{x) xo{x) , (B.4) 



30 
25 
20 
15 
10 
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2 . 5 



6x 



Figure Al. Percentage difference between the numerical coefficients of ei £2 
in P/j from WKB* and GFM as a function of 6x- The difference vanishes for 
6t — 1.93 (~ 2, the value used in the main text). 
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and using the standard formulae for second order differential equations with a source 
term, so that 



9n{x) 



xliy) 



(B.5) 



The problem has therefore been reduced to computing 2 n integrals. 

The above formalism can be modified in order to provide an alternative way for 
calculating the power spectra, spectral indices and runnings to second order in the 
HFF. We start from Eq. © and ifTUI) and the exact relations 



■ £2 



zs \ ' 2 ^ 2 ^ ^ 4 

With a redefinitions of the wave-function and variable 



1 



1 



■ £2 £3 



we obtain 



where 



X = In (— fc 77) 

V^i) = , 



(B.6) 

(B.7) 
(B.8) 



6'"= - ^-7]'a'H' (2-ei) 



■ £2 



■ £1 £2 



On now using the approximate relation 

V- 77 (l+€l+€?+€l£2) 

an 



the frequencies become 



T + 3£i 



4e 



4 6162 



£2 £3 



:e-- (- + 361 + -62+46? + _e,e2 + _e2 



/2x 



2£2£3 



(B.9) 



(B.IO) 



(B.ll) 



where we only keep the second order terms in the (still i-dependent) HFF. We can 
now expand linearly each HFF around the horizon crossing to obtain 



+ 



d 



e — 



e — 



9 o 

- +3ei +46^ + 46162 



V'T = —3 61 62 X IpT 



(B.12) 



9 o ^ 3 . ^ .2 13 . . 1 .2 1 ^ ^ 

-+3 6i + - 62 +4 6i + y 6162 + 4 ^2 + 3 ^2£3 
1 



-3 ( 61 62 + - 62 63 ) , 
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which, beside the use of St instead of x and the appearance of uj instead of w, is 
equivalent to the expansion in Eq. (|B.1|) for / = (D^ and can be solved by means of 
a standard perturbative approach (such as in Ref. [S])- The method is close to GFM 
and gives the same results of Eq. (41) in Ref. [3] and Eq. (29) in Ref. [201 for the scalar 
and tensor power spectra, respectively. 

Let us finally note that, for the GFM, the scalar spectral index and its running 
are calculated by taking the derivative of Eq. (41) in Ref. i.e. the power spectra 
evaluated aX k = aH, in terms of Infc, while in our MCE (but also in the WKB 
approximation 'S', '^) the dependences in ln(fc/A:») are obtained explicitly (see, for 
example, Eqs. H32a|l and (j32b|) '). 
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